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Variability of a stellar corona on a time scale of days 
Evidence for abundance fractionation in an emerging coronal active region 

R. E. Behar^ and S. A. Drake"* 



(N 
O 

o 

Q: 

(N 



O 



> 

O 



X 



Max-Planck-Institut fiir Extraterrestrische Physik (MPE), Postfach 1312 85741 Garching Germany 

School of Physics and Astronomy, The Raymond and Beverly Sackler Faculty of Exact Sciences, Tel-Aviv University, Tel-Aviv 
69978, Israel, e-mail: nordon@astro . tau . ac. i 1 

Physics Department, Technion Israel Institute of Technology, Haifa 32000, Israel, e-mail: beharOphysics . technion. ac . il 
NASA/GSFC, Greenbelt, Maryland 20771, USA. e-mail: Stephen . a . drakeOnasa . gov 



Received...; accepted.. 



ABSTRACT 



Elemental abundance effects in active coronae have eluded our understanding for almost three decades, since the discovery of the First 
Ionization Potential (FIP) effect on the sun. The goal of this paper is to monitor the same coronal structures over a time interval of 
six days and resolve active regions on a stellar corona through rotational modulation. We report on four iso-phase X-ray spectroscopic 
observations of the RS C Vn binary EI Eri with XMM-Newton, carried out approximately every two days, to match the rotation period 
of EI Eri. We present an analysis of the thermal and chemical structure of the EI Eri corona as it evolves over the six days. Although 
the corona is rather steady in its temperature distribution, the emission measure and FIP bias both vary and seem to be correlated. An 
active region, predating the beginning of the campaign, repeatedly enters into our view at the same phase as it rotates from beyond 
the stellar limb. As a result, the abundances tend slightly, but consistently, to increase for high FIP elements (an inverse FIP effect) 
with phase. We estimate the abundance increase of high FIP elements in the active region to be of ~75% over the coronal mean. 
This observed fractionation of elements in an active region on time scales of days provides circumstantial clues regarding the element 
enrichment mechanism of non-flaring stellar coronae. 

Key words. stars:atmospheres - stars:coronae - stars: abundances - stars:individual:EI Eridani - X-rays:stars 



1. Introduction 

The best studied cool-star corona is that of our own sun, which 
also provides the only spatially resolved stellar coronal source 
of Extreme Ultraviolet (EUV) and X-rays. From X-ray images 
of the sun we know that the solar corona in its active state is 
non-uniform and that a significant fraction of the emission in 
these bands comes from magnetic loop structures concentrated 
in active regions around photospheric sun-spots. Evidently, these 
structures evolve in brightness on time scales of hours. The mor- 
phology as well as the metal abundances vary on longer time 
scales of days. Coronal structures on the Sun may survive for 
a few weeks, although tracking them continuously from the 
ground is difficult due to the ~ 24.5 day solar rotation period. 

Active coronae of other cool stars can be up to five orders of 
magnitude more luminous and an order of magnitude hotter than 
the solar corona. Although unresolved, they show many charac- 
teristics similar to those of the solar corona, and are therefore 
believed to be structured and to evolve in a similar way. For the 
most part, late-type stars have envelopes in which convection 
and differential rotation interact to drive a dynamo, which pro- 
duces magnetic flux tubes. The flux tubes float into the upper 
atmosphere and eventually the corona. It is generally unclear, 
however, whether stellar coronae actually have larger covering 
of their stellar photosphere, or even qualitatively different mor- 
phologies, such as inter-binary magnetic structures. Stellar coro- 
nae show high variability, often attributed to continuous micro- 
flaring, which is observed on the sun, but perhaps with different 
statistics (e.g.. lCaramazza et al ] l2007h . 



Solar coronal chemical abundances are different from the 
solar photospheric composition. The solar corona itself is non- 
uniform in that respect. While in coronal holes, the abundances 
are close to photospheric, abundances in active regions and 
in the solar wind are significantly different. The abundance 
p atterns depend on t he first i onization potential (FIP) (Meve^ 
119851 iFeldrnanli 19921; ISheelevllT995h . Elements with FIP<10 eV 
(low FIP) are enriched in the corona relative to elements with 
FIP>12 eV (high FIP). Averag ing over the entire solar disk, 
the enrichment factor is about 4 f Feldman & Widirig||2002l and 
references therein). Such a FIP bias, in varying intensities, has 
been reported for many coronal structures, from closed mag- 
neti c loops to streamers (see reviews by iKohl et all l2006t 
Dos chek & Feldmanll2010h . Currently there is no generally ac- 
cepted model that explains the FIP bias. 

The FIP effect became even more puzzling when high res- 
olution X-ray spectra from XMM-Newton and Chandra revealed 
that active stellar coronae do not follow the solar FIP pattern. In 
some cases, the high-FIP elements are even enriched compared 
to th e low-FIP ones, an eff ect labeled the inverse-FIP (IFIP) ef- 
fect dBrinkman et al.ll200 1'). Later studies revealed that FIP and 
IFIP biases are correlated with coronal activity (and age): Highly 
active stars sh ow an IFIP effect, while less active coronae have a 
solar FIP bias jAudard et al. 2003; Te lleschi et al. 2005) . imply- 
ing a transition on stellar-evolution time scales of Gyrs. 

On the face of it, the correlation with activity supported 
previous suggestions tha t abundances may be linked to flares 
(e.g., Gudel et al. 1999i lOsten et all I200I iGudel et al.1 l2004t 
iNordon et al.u2006i) . However, direct systematic measurements 
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of abundance variations during flares show that regardless of 
coronal FIP/IFIP bias, during large flares the abundances tend to 
photospheric values (Nordon & Behar 2007, 2008). The mecha- 
nism for FIP/IFIP is, therefore, related either to the long-term 
evolution of coronal structures, or to the continuous micro- 
flaring that produces the quiescence emission, but not to large 
flares. 

In recent years, a model has been proposed in which abun- 
dance fractionation is driven by m agneto hydrodynamics (MHD) 
waves ( lLamind2004 120091 120 1 21) . This 'pondermotive' force is 
controlled by the energy density of M HD waves at the bases of 
the coronal loops. IWood et alJ (l2012h studied a sample of X-ray 
faint (Lx < 10^^ erg s"') dwarfs and found no correlation with 
X-ray luminosity or surface flux, but a good correlation with 
spectral type instead. According to this relation M dw arfs are 
IFIP biased while early G dwarfs are FIP biased. Wo od et aTl 
(l20T2h claim that this shift from FIP to IFIP can be explained 
by the inversion of the pondermotive force at certain conditions. 
The latter study is limited to X-ray faint dwarfs and may not be 
applicable to the typical highly active binaries from other stud- 
ies that are more than two order of magnitude more luminous in 
X-rays, emitted mostly by a giant or sub-giant star (RSCVn type 
systems). 

Despite the assumed similarity between solar and stellar 
coronae, we have lacked a study of the evolution of a stellar 
corona on time scales of days, which is the relevant time for 
coronal loop evolution on the sun. The goal of the observations 
presented in this paper is to seek and measure the variability of 
temperatures and abundances in an active corona over several 
days. In order to avoid confusion between rotational variability, 
which is inevitable due to the rapid rotation typical of active stel- 
lar targets and the corona being composed of discrete active re- 
gions as on the sun, and genuine coronal-activity variability, the 
present observing campaign is designed to return to the target 
each time at exactly the same orbital phase. 

1.1. ElEri 

In order to carry out a few iso-phase observations, we chose 
flie target of EI Eridani (EI Eri, HD 26337), an X-ray bright 
(Lx = 10^''^ erg/s) RS CVn binary, whose orbital period is 1.95 
days (168.5 ks), which is very similar to the exactly 2.0 day or- 
bital period of XMM-Newton around Earth. This enabled XMM- 
Newton to observe the exact same orbital phase of EI Eri dur- 
ing four consecutive orbits, and for sufficiently l ong exposures 
each time (~ 35 ks). EI Eri is a nearby (50.8 pc, Ivan LeeuwenI 
l2007h . rapidly rotating (vsin/ - 51kms '), non-eclipsing, rota- 
tionally locked, close RS CVn bi nary with an inc lination angle to 
the rotation axis of ~ 56° + 4.5° (IWashuettl et al. 2009a). These 
properties make EI Eri a popular target for Doppler imaging and 
starspot monitoring that indicate the primary star to be highly 
active. EI Eri consists of a solar-mass G5IV primary and a much 
fainter M4-5 secondary, with a possible 19-year orbit of a third 
body (IWashuettl et al.ll2009ah . 

Long term photometric monitoring of EI Eri seemed to 
indicate a ~11 yr brightness cycle that was conjectured to 
be an analog of the solar activity cycle (Olah & Strassmeier 
I2002h. but unlike the sun, l ong term Doppler imaging of EI Eri 
by IWashuettl et afl (l2009bl) shows little va riability of i t s cool 
starspots over 1 1 yrs. As for the hot corona, lOsten et alj (|2002|) 
observed EI Eri with ASCA (and EUVE) and fitted its spectrum 
with a 2T model of kT - 0.72 keV and 2.2 keV, but the abun- 
dances could not be tightly constrained with ASCA. The XMM- 
Newton CCD spectra of a short flare on EI Eri on 2009-08-06 



during the present campaign was analyzed by iPandev & Singh! 
(1201 2l) . who suggested the flare originated from a loop the size 
of ~ Rt (primary), which is about a tenth of the binary separation 
distance. 

2. Observations, Data, and Analysis 

In Table 1 we list the log of the four XMM-Newton observations 
of EI Eri. Note the start times for each observation, which are 
separated by approximately the EI Eri orbital period of 168.5 ks. 

2.1. Data Reduction 

XMM-Newton data were reduced using the Science Analysis 
Software (SAS - ver. 10.0.0) in the standard processing chains 
as described in the ABC Guide to XMM-Newton Data Analysifl 
and using the most updated calibration files (date 22.06.2010). 
Pileup was checked using the epatplot routine of the SAS soft- 
ware package. No significant pileup was detected in either of the 
EPIC instruments and thus events from single to quadruple were 
retained. Background was low during all the observations and no 
significant amount of time was lost. 

2.2. Spectral Analysis 

We measure the temperature and abundance properties of the 
EI Eri corona by fitting thermal plasma models to the observed 
spectra. XMM-Newton comprises several different instruments. 
The longer wavelengths (6-38 A) are covered by the two high 
resolution reflection grating spectrometers (RGS) that resolve 
the emission lines. The high energy end up to 15 keV is covered 
by the two EPIC-MOS and single EPIC-pn CCDs, that provide 
lower spectral resolution. 

In order to study the thermal and chemical composition of 
the source plasma, it is imperative to first properly account for 
its Emission Measure Distribution that is defined as 

EMD{T) = nenudVldT. (1) 

where Ug and uh ait the electron and proton number densities, 
respectively, and dV/dT is the temperature distribution of the 
plasma volume. The EMD is usually characterized by several 
discrete temperature bins. Here, we use seven temperature bins 
with identical abundances. To allow for a uniform analysis of 
the four observations, the temperatures of the middle five com- 
ponents kT = 0.7, 1.0, 1.5, 2.0, and 2.5 keV are held fixed, 
while the lowest and highest temperatures are allowed to vary 
and are determined by fitting the model to the observed spec- 
tra. We fit the observed spectra using the XSPEC ("A rnaudl 19961) 
software package and a VAPEC ( Smith et al. 2001 ) plasma spec- 
tral model. For high statistics, we fit the two RGS, EPIC-pn, and 
EPIC-MOS 1 spectra of each observation (orbit) simultaneously. 
EPIC-M0S2 was used in timing mode and its spectra are sys- 
tematically offset by a small amount from EPIC-pn and EPIC- 
MOS 1 that agree with each other very well. For this reason, we 
exclude it from the spectral fitting. We ignore the very low en- 
ergy (E < 0.35 keV) portions of the CCD spectra. Visual inspec- 
tion of each RGS spectrum indicates the models do an overall 
good job of fitting the lines. See the appendix for the spectra and 
fitted model. The absorption colum n density towards EI E ri is 
negUgible (Nh = 1.2 x 10^" cm ^ iPandev & Singhll2012l) and 
was not included in the model. 

' See http://heasarc.gsfc.nas a.gov/docs/xrmn/abcl 
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Table 1. Observing log 

Obs. ID Start time End time Exposure [ks] 

060317201 2009-08-01X00:02:50 2009-08-01X09:13:05 32.84 

060317301 2009-08-02X23:06:16 2009-08-03X09:24:46 36.92 

060317401 2009-08-04X22:20:49 2009-08-05X07:42:41 33.54 

060317501 2009-08-06X20:59:39 2009-08-07X07:31:30 37.48 




0.00 0.05 0.10 0.15 0.20 

Phase 



Fig. 1. EPIC-pn light curves during the four consecutive EI Eri 
rotation orbits. The time axis has been folded to represent the 
rotation phase within the 1.95 days period of EI Eri. 



The seven plasma components provide only a general idea of 
the temp erature distributio n in the plasma. It is known (see most 
recently, iTesta et £0120 12h that the EMD reconstruction should 
be treated with caution. More specifically, the EMD is highly 
degenerate between adjacent components. Therefore, we keep 
track of the covariance matrices (between parameters), which are 
fully taken into account when propagating the errors of all de- 
rived quantities. Due to this degeneracy, the more important and 
reliable quantity is the cumulative, or progressively Integrated 
Emission Measure which we denote by lEM. The 

calculated errors for the lEM are more meaningful than those 
of the EMD, as they are essentially independent of the choice 
of temperature bins, as opposed to errors on the EMD that are 
affected by the arbitrary choice of temperatures and by the prox- 
imity of the other temperature components. Monte Carlo (MC) 
simulations are used for error estimates. The MC simulation ran- 
domly rolls a new set of EMD values according to the correla- 
tion matrix and individual-component errors obtained from the 
best- fit solution. The central 68% of the distribution for each 
point in the lEM is taken as the equivalent of the 1-cr confidence 
interval. For a more detailed description of the error propagation 
see Appendix|B] 

3. Results 

3.1. Light Curves 

By synchronizing the XMM-Newton exposures with the rotation 
phase of EI Eri, we are able to observe the same stellar surface 
during each visit, and to thus remove most of the rotation ef- 
fects. The EPIC-pn light curves of the four exposures are plotted 
in Figure [U versus the rotaional phase. Only about 20% of the 



orbit, i.e., 35 ks or ~0.4 days, is covered during each visit, and 
consecutive visits are separated by the remainder of the EI Eri 
period (~2 days). 

Between phases - 0. 1, all 4 exposures show a steady count 
rate with some variability on time scales of minutes, but at a 
low amplitude that is not much above the Poisson noise. This 
variability is typical of active coronae, and has been ascribed to 
micro-flaring. The second half of each observation, i.e., phases 
0.1 - 0.2, features a gradual rise in flux that is interpreted as the 
recurrent appearance of an active region on that part of the star 

This brightening with phase is most dramatic during the third 
observation (401), ~4 days into the monitoring, but should still 
not be confused with a coronal flare that would tend to sharply 
rise within a few ks and quickly decay over ~ 10 ks. Indeed, 
during the last exposure (obs. 501), a large flare is superimposed 
on the gradually increasing light curve, rising sharply (~1.5 ks) 
about half way through the exposure and decaying exponentially 
on a time scale of ~4.8 ks. For analysis purposes, we therefore 
further split observation 501 into its quiescence and flare peri- 
ods, as indicated in Fig. [1] and refer to them as obs. 501Q and 
50 IF, respectively. 

The mean count rate between visits (~ 130 ks) also varies, 
increasing within the first 4 days, between obs. 201 and 401 by 
~ 20%, and then decreasing by ~30% by the beginning of obs. 
501 . Interestingly, just before the flare and ~5 days into the mon- 
itoring, the steady quiescent level is observed to be back to its 
initial level when the monitoring started (obs. 201, 301). 

3.2. Temperature 

In Figure|2]we plot the EMD(T) for each observation based on 
our 7-temperatures best-fitted models, where each EMD bin is 
due to a single temperature component. The width of each plot- 
ted bin is arbitrary and designed to give an impression of a con- 
tinuous distribution, although in reality these are iso-T compo- 
nents. As with the light curves, the EMD of observations 201, 
301, and 501Q are fairly similar, primarily showing emission of 
plasma up to kT ^ 2 keV, while the EMD of observation 401 
features some excess high-T emission up to kT ^ 4 keV. During 
the flare (obs. 501F), significant emission measure is observed 
at kT > 4 keV. 

The lEMs for each observation are then plotted in Fig. [3] 
The lEM in observations 201, 301 and 501Q are consistent with 
no variation in the total emission measure (see the high tem- 
perature end of Fig. [3]), and show very small variations below 
1 keV. On the other hand, the EI Eri corona during observation 
401 has less emission measure below 1 keV, but a higher total, 
implying a significant excess of plasma at higher temperatures, 
mostly between 2.5-3.5 keV. During the short flare (obs. 501F), 
one can see an excess of emission measure between 1-2 keV and 
another small component above 3 ke V that brings the total emis- 
sion measure to a level ~ 40% higher than the quiescent states of 
obs. 201, 301, and the pre-flare 501Q, which is also somewhat 
lower than during obs. 401. In short, the emission measure of 
the quiescent corona of EI Eri can vary between rotation orbits. 
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Fig. 4. Variability of the EMD-weighted mean temperature. 
Horizontal error bars indicate the exposure time of each obser- 
vation. 



Fig. 2. Emission measure distribution as obtained by the best- 
fitted 7T models of each observation. 
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Fig. 3. Cumulatively integrated emission measure distribution 
IBM calculated from the EMDs of Fig. |2] Errors take into ac- 
count the degeneracy between the different temperature compo- 
nents of the models. 



and on time scale of days, but mostly in the 1 -2 keV tempera- 
ture regime. Excess emission measure of hotter plasma up to 4 
- 5 keV is observed both in active regions (obs. 401), and more 
dramatically during the flare (obs. 501F). 

A concise way to quantify the plasma conditions is through 
its mean temperature < kT >, weighted by the EMD of each 
component. The errors are again computed using an MC simula- 
tion, while taking into account the correlations between the fitted 
EMDs of the various model components. The resulting < kT > 
values are plotted in Figure |4] for each observation and against 
time elapsed from the beginning of the first observation. A sim- 
ilar pattern as before emerges. When the source is brighter (e.g., 
Obs. 401) it is also hotter on average. Obs. 201 and 301 do not 
strictly follow this pattern, but the temperature decrease from 
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Fig. 5. Coronal abundances relative to Fe and expressed relative 
to the solar X/Fe abundance ratio jAsplund et al.ll2009t) . Points 
from different observations are progressively shifted to the right 
along the horizontal axis in the order in which they were taken 
(~2 days apart) to enhance visibility and to emphasize the time 
variability. 



Obs. 301 to Obs. 201 is insignificant. In terms of < kT >, we 
observe temperature variations on time scales of days, and most 
dramatically before and during the flare. Temperature variability 
during quiescence is more subtle, with a mean kT of 1 .60 keV 
and a standard deviation of 0. 13 keV. 

3.3. Abundances 

In Figure|5]we plot the relative abundances of the major elements 
with respect to Fe, in solar-abundance units dAsplund et al.l 
'2009!), and as a function of their first ionization potential (FIP). 
In the plot, the data points are progressively shifted along the 
X-axis according to the order in time in which they were mea- 
sured, thus creating a time-like sequence for each element. The 
abundances are expressed relative to the Fe abundance instead 
of the more common absolute values relative to H. The absolute 
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Fig. 6. FIP bias (FB from Eq. |2]i for each observation and as 
a function of time. Horizontal error bars indicate the exposure 
time of each observation. 



Fe abundance with respect to H is 0.32, 0.32, 0.38 & 0.31 solar 
in observations 201, 301, 401 & 501Q respectively, with a typ- 
ical error of +0.03. In the APEC model and the global spectral 
fit approach in the XSPEC software that we use, the absolute 
abundance is constrained by the flux in the Fe lines and the over- 
all fit to the continuum. While lines are emitted from a limited 
range of temperautes, continuum flux at all RGS wavelengths 
is contributed from all temperature components. Absolute abun- 
dance is therefore sensitive to the underlying temperature distri- 
bution, which is not so well constrained at the high temperatures. 
Conversely, estimates of abundances relative to Fe rely on well- 
measured line ratios that tend to emerge from a limited and over- 
lapping range of temperatures. The large number of observed Fe 
ions and lines, whose peak emissivities span a wide tempera- 
ture range from kT ~ 0.3 - 1 .5 keV makes these measurements 
significantly less sensitive to the temperature distribution in the 
plasma. The selection of the reference element for describing 
abundances does not affect the general pattern of abundance vari- 
ations. We normalize the abu ndance ratios t o those of the solar 
photospheric composition of ' Asplund et alj (l2009h in order to 
plot them on the same scale and to discuss the FIP effect. 

As seen in Fig. |5] EI Eri shows the t ypical inverse F'lP ef- 
fect observed in the most active coronae ( Brinkman et al. I I2OOII 
lAudard et al.ll2003h . In particular, the Ne abundance is enriched 
by a factor ~3 compared to Fe, and a lesser enrichment is ob- 
served for C, N, and O (all high FIP) over Fe (low FIP), The 
abundance ratio of the low-FIP elements Mg and Si to Fe is ap- 
proximately solar, although with respect to H, they are all sub- 
solar (~ 0.3). 

Little variation of the abundances with the level of activity 
is observed. The O (high FIP) abundance does seem to be low 
during the high-states of obs. 401 and 50 IF (flare), and high dur- 
ing the low-activity state of obs. 501Q. This is reminiscent of the 
chromospheric evaporation effect of solar-composition material 
during flares in an oth erwise (quiescent) inverse-F IP coronae, 
which was observed bv lNordon & Behad (l2007ll200 8> m a sam- 
ple of flares. The highly enriched Ne abundance is slightly re- 
duced in the high state of obs. 401. However, despite the good 
signal, no significant variation in the Ne abundance can be de- 
tected during any other period, not even during the flare, which 
would be expected if chromospheric evaporation takes place. 



The Mg/Fe abundance ratio comes down by ~ 20% during the 
flare, but this can not be a FIP related effect, as both elements 
are low-FIP. 

In order to better quantify the abundance variations, a single 
nu mber that quantifies the FIP bias can be defined (following 
No rdon & Behadl2008l) 

FB = <log(Az/Aze)FIP<10eV> - <log(Az/Azg)FIP>10eV> (2) 

where {log(Az/Azg)) is the error- weighted mean abundance (in 
log) of the low/high FIP elements, relative to a set of reference 
abundances, in this case the solar abundances. We use the tradi- 
tional solar distinction between low FIP (<10 eV) and high FIP 
(>10 eV). A positive FB value indicates a solar-like FIP bias, 
while a negative FB value indicates an IFIP bias. 

The FIP biases in the five observations (20 1 , 30 1 , 40 1 , 50 1 Q, 
501F) are shown in Fig. |6] The FB variations over the course 
of the present campaign are not dramatic, but still some trends 
can be identified. The FB in Fig. |6] generally varies with activ- 
ity, namely the FB increases, i.e., tends to (solar) photospheric 
values, as the activity increases, rising from obs. 201 to obs. 
301, and further in obs. 401, and then decreasing in obs. 501Q 
when the count rate returns to its lower level. For comparison, 
see the light curves in FiglU This trend breaks down during the 
flare when the activity strongly increases, but the FB decreases. 
At face value, this means that as the (quiescent) coronal ac- 
tivity increases, perhaps more small active regions appear, and 
more high-FIP elements are preferentially brought up into the 
corona. This process is reminiscent of the evaporation of chro- 
mospheric material observed in the largest flares on the most 
active (strongest IFIP) coronae (Nordon & Behar 2007, 2008)). 
On the other hand, the present flare on EI Eri does not follow 
this trend, but the weak abundance variation observed is consis- 
tent with the abundance variations observed in flares on coronae 
with similar quiescent FB values of about -0.2 (see Fig. 5 in 
Nordon & Behar..2008.) . 

4. Rotation Effects 

The four exposures that all cover the same rotation phase al- 
low us to study the dependence of temperatures and abundances 
on phase, by further breaking up each observation into smaller 
segments and co-adding the similar phase segments from each 
observation. If active coronal regions are compact, and sparse 
in longitude, then the rotation would create variations in the ob- 
served X-ray emission properties, as active regions rotate in and 
out of view, over the limb of the star. If the active regions last 
longer than the monitoring time of 6 days (as solar active regions 
do), then stacking data from the same phase, but from different 
orbits, can increase the signal of this rotational modulation. 

During each full observation of ~36 ks (Table 1), EI Eri ro- 
tates by 72 degrees, i.e., 20% of the orbit. The inclination an- 
gle of the rotation axis to the line of sight of 56° implies only 
about 10% of its surface gradually appears during the observa- 
tion and a different 10% gradually moves out of sight. We di- 
vide each observation into four quarters of ~ 9 ks each, during 
each of which the EI Eri system rotates by 18 degrees. The data 
are then reduced in quarters and the spectra of each instrument 
are combined per phase bin. The two late phase bins of Obs. 
501 are excluded from the analysis, due to the flare (see Fig.lTJ. 
Subsequently, we perform the same EM and abundance analysis 
as before, but for each phase-stacked spectra. 

Figure Q shows the IBM for each phase bin, Ql - Q4. The 
first three phase bins show the same temperature structure. In 
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Fig. 7. Cumulatively integrated emission measure distribution Fig. 8. Mean temperature for the four rotation phase bins. Very 



(I EM) as a function of the four phase bins Ql - Q4. 



little evolution is observed compared to the days-scale tempera- 
ture variations (Figure|4]i. 



the last phase bin (Q4), there is an excess of emission at temper- 
atures below ~1.5 keV, which slightly increases the total lEM 
of Q4, although it is just consistent with the previous quarters 
to within the errors. This excess emission measure is a result of 
the rise in flux observed in the light curves (Figure[T]) during the 
last quarter of all of the observations, but most conspicuously 
in Obs. 401. A similar rise in flux can be seen also in the last 
quarter of Obs. 501, which is somewhat masked by the flare, 
and therefore excluded from the analysis. Figure |7] can be ex- 
plained by an active region that is not visible at all during the 
two earlier quarters, that comes into view from beyond the limb 
during Q3 or Q4 (i.e., after ~ 40 degrees of rotation), and that 
gradually brightens during the four days after the beginning of 
the XMM-Newton monitoring, i.e., it becomes particularly bright 
during Obs. 401 and Obs. 501. 

In Figure [8] we plot the mean temperature as a function of 
phase. The value of kT^g is very stable between the observed 
rotation phases, despite the added emission measure. In fact, it 
shows much less scatter than the scatter from one observation 
to the other over several days (Figure |4|. The lack of temepra- 
ture variation is somewhat expected, since each phase represent 
a time-average of 4-6 days, and variations on shorter time scales 
tend to average out. Note that here too, the two last quarters of 
Obs. 501 are excluded from the analysis due to the flare. The 
emerging active region (presumably a magnetic loop) reaches 
the mean temperature of (kT ~ 1 .5 keV), which is very similar 
to the mean temperature of the coronal regions on the rest of the 
EI Eri observed area. 

The mean abundances for each rotation phase bin are plotted 
in Figure|9] Again, the data points are progressively shifted along 
the horizontal axis, in order to give an idea of the evolution with 
phase. The variation in abundance with phase for each element 
is very small, and statistically insignificant between consecutive 
bins. However, Mg, O and Ne, the three elements with the best 
signal do appear to vary monotonously with phase. The abun- 
dance ratio of Mg/Fe decreases while that of O/Fe and Ne/Fe 
increases. The other elements tend to show a similar monotonic 
variation, although with much larger uncertainties. While the 
trend is of low signal in each element, the repeated pattern in 
several elements makes it more believable, and less likely due to 
noise. Moreover, since the trend of the low-FIP Mg is opposite 



to that of the high-FIP Ne and O, we suspect it is due to a FIP 
related efifect and not due to only Fe changing with phase. 

As in Sec. 13.3! we wish to enhance the individual abundance 
signals as to reflect potential FIP related effects by computing 
the FIP bias from Eq. |2]for each phase-stacked spectra. The FIP 
bias as a function of phase is plotted in Fig.[TOl The FB shows 
a mild, yet monotonous decrease with phase. This can be inter- 
preted as the emerging active region contributing an inverse FIP 
bias. In other words, the active region (Q3, Q4) has slightly more 
high-FIP elements than the background corona (Ql, Q2). Note 
that this time-averaged rotation effect is different from the time- 
evolution observed over days (Fig. 6), in which the FB increases 
with increased activity. Here, the increased activity coming from 
a specific region on the star is the likely cause of the decrease in 
FB. 

The weakness of the observed signal is not unexpected, as 
the fraction of the EI Eri surface rotating in and out of our view 
is not very large, ~ 5% between quarters. Say we measure dur- 
ing rotation phase Q4, in which the active region is visible, a 
mean abundance A^^ which deviates from the mean abundance 
during phase Ql A^' (in which the active region is not visible) 

by ff « 0.1, i.e. A^^ = (1 + a)A^^\ The total count rate is in- 
creasing with phase, indicating that this efifect is likely due to 
a different abundance in the active region which is rotating into 
view, for which A^*'"" = (1 -H I3)A^2 ■ Since the contributions to 
the spectrum of both the background and the active-region abun- 
dances are weighted by their respective emission measures, and 
we attribute all the excess emission measure during Q4 over Ql 
to the active region, then one can easily show that; 



^ 1 -(£M(i)/£M(4)) 

where EM*'*) and EM*" are the total emission measures in Q4 
and Ql respectively. According to Figure |7] the EM ratio is 
EM^'*''/EM^ X 1.15. Therefore, the actual abundance variation 
in the active region by comparison with the spatial-mean back- 
ground corona is at the /3 « 0.75 level, or a 75% stronger IFIP ef- 
fect than the mean corona. Admittedly, this conclusion for EI Eri 
is based on a marginal signal and is, thus, associated with ap- 
preciable uncertainty, but it does point to significant abundance 
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Fig. 9. Coronal abundances relative to Fe and expressed relative 
to the solar X/Fe abundance ratio ("Asplu nd et al.ll2009h . Points 
from different phase bins Ql - Q4 are progressively shifted to 
the right along the horizontal axis to enhance visibility and to 
emphasize the variability with phase. 




0.10 0.15 

Rotation Phase 

Fig. 10. FIP bias (FB) as a function of rotation phase aver- 
aged over all available observations in four iso-phase intervals. 
Horizontal error bars reflect the ~9 ks duration in each bin. 



variations in a compact active region with respect to the mean 
corona. 

5. Discussion and Conclusions 

The light curves and abundance patterns in the corona of EI Eri 
show an active region appearing over the limb as the corona 
rotates with the star. Although on time scales of days, the in- 
crease of the activity level seems to lead to an increase of the FB 
(i.e., less IFIP bias, abundances more akin to the photosphere), 
the appearance of the active region with phase features an in- 
crease in high-FIP elements and a lowered FB (i.e. a tendency 
to a stronger IFIP bias). Two different effects are at play: one is 
the evolution of the emission measure and coronal abundances 
on a time scale of days between exposures, spatially integrated 
and averaged over the entire visible surface of the star The other 
is an already-existing, compact and bright active-region that per- 



sists over the six days duration of the monitoring campaign and 
rotates into view towards the end of each exposure. Thus, the 
active regions is manifested as variations with phase instead of 
time. 

One cannot tell whether the days-scale variations in count 
rates, emission measure, and abundances are due to the evo- 
lution of existing coronal loops, or due to the appearance of 
new coronal loops and disappearance of old loops. However, 
the correlation in which an increase in emission measure and 
count rates corresponds to a decrease in the IFIP bias (increase 
in FB, closer to photospheric composition) offers an interesting 
interpretation. If the brightening is due to the emergence of new 
loops, and those appear initially with nearly photospheric com- 
position, then the mean coronal abundances will tend towards 
photospheric composition during these times. Conversely, dur- 
ing times in which no new loops emerge and existing loops 
age, evolve in composition and perhaps some fade, the mean 
abundances will tend to resemble the composition of the older, 
evolved loops, i.e. different than photospheric. 

Indeed, on the sun, it has been reported that solar coronal 
loops evolve on time scales of days (e.g., Widing & Feldma^ 
120011) . Newly emerging structures tend to have photospheric 
composition, but develop the solar FIP bias within a few days. 
Moreover, on the sun, the typical enhancement of low FIP abun- 
dances relative to photospheric values by a factor of 4 is reached 
within 2-3 days and continues to evolve further, reaching a fac- 
tor of 8 or more. The analogous interpretation and one in agree- 
ment with the present observations is that new coronal loops on 
EI Eri emerge with photospheric composition and develop an 
IFIP bias over time scales of days. The spatial mean IFIP bias is 
determined by mean age of the loops. Thus, periods of gradual 
brightening of the corona are characterized by the emergence of 
young loops that reduce the global IFIP bias, while periods of 
gradual dimming are characterized by aging of existing loops 
that increase the global IFIP bias. Note that the negative correla- 
tion between emission measure and IFIP bias as described above 
relates to the days scale evolution of a given stellar corona and is 
not at odds with the general positive correlation found between 
activity level and tendency to IFIP bias in stellar populations. 

Rotational mod ulation in X-rays and EUV has been re- 
ported for HR 1099 jAgrawal & VaidvalI988tlDrake etalFm^ 
lAudard et al.l l200lh . from light curves, but phase dependent 
spectroscopy was not possible, due to the constant flaring of that 
source, which makes it difficult to disentangle between rotational 
modulation of active regions and their intrinsic activity variabil- 
ity. A different, yet notabl e, phase depe ndent measurement is 
that of line profiles by Bric khouse et al.l (j2001), which enabled 
latitude determination of active coronal regions on the stellar bi- 
nary 44i Boo. The present measurement, however, is the first to 
measure phase dependent abundances. 

The active region that appears towards the end of each ex- 
posure persists throughout the six days of the observing cam- 
paign. It therefore likely contains loops that are at least a few 
days old. In the scenario which we proposed above, these loops 
will have already evolved (over days, possibly starting even be- 
fore the beginning of the campaign) and their abundances devel- 
oped a strong IFIP bias, stronger than the coronal mean bias that 
includes younger loops. Thus, as this region comes into view, 
the spatially averaged abundances tend towards a stronger IFIP 
bias (Figure [Toll. 

After showing that large flares tend to counter the coronal 
FIP eff'ect (Nordon & Behar 2007, 2008) - the present work pro- 
motes the notion that the abundances in individual coronal ac- 
tive regions, differ from the mean coronal abundances as well. 
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It is not clear whether large flares show coronal abundances that 
are closer to photospheric composition than the quiescent corona 
because of added material from the photosphere/chromosphere 
(chromospheric evaporation) or because these flares happen 
preferentially in young active regions where the abundances are 
closer to photospheric than the global mean. The abundances 
measured during the flare observed in this campaign show an in- 
crease IFIP bias (more negative FB, Figure|6]l instead of the typ- 
ical tendency to bring abundances towards photospheric com- 
position. If the flare is related to the limb active region that is 
strongly IFIP biased, it is possible that the mean abundances, 
weighted between the flaring active region, evaporated chromo- 
spheric material and surrounding quiescent corona end up being 
dominated by the abundances of the active region. 

While the proposed scenario in which young loops emerge 
with photospheric-like composition and then develop over time- 
scale of days to an IFIP bias can explain our observations, it does 
not explain why the abundances evolve an IFIP eff'ect instead of 
a solar-li ke FIP effect. For a theory t hat attempts to resolve this 
issue see iLamingl (|2004 l2009l |20 1 2h . 
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Appendix A: Spectra 

Sample spectra and best-fitted model are presented in this sec- 
tion. The X-ray spectra obtained from our first observation ID 
060317201, as measured by all XMM-Newton instruments used 
in the analysis are shown in Fig. I A. II This observation has the 
lowest count rates (see Fig. [TJ. Consequently, the spectra of 
the three other observations have slightly better signal-to-noise 
ratios, but are overall similar to those shown in Fig. lA.ll In 
Section |4] each observation is divided into four segments, and 
the four iso-phase segments are then stacked, resulting in spec- 
tra that are, again, with numbers of photons and signal-to-noise 
ratios that are comparable to those of the individual observations. 

Appendix B: lEM and <kT > Error Estimates 

The EMD solutions to coronal emission tend to be highly degen- 
erate as emission measure can be traded between different, but 
close temperatures with little effect on the resulting spectrum. 
Thus, when representing the EMD with isothermal components, 
the uncertainties in the normalizations of the various compo- 
nents are highly coiTelated, especially when many components 
are used. The correlation tends to be negative and increases as 
the components are closer in temperature, i.e., adjacent temper- 
atures can compensate for each other. This results in large errors 
on each individual temperature component, which can be seen 
in Fig. 12] On the other hand, the exchange of emission measure 
between close-temperature components has only a small effect 
on the cumulative, or integrated emission measure {IBM ). Still, 
the spectral line emissivities are temperature dependent and thus 
the IBM is not perfectly conserved. In addition, there are the 
uncorrected uncertainties associated with the limited statistics 
of the data. In order to estimate the highly nonlinear contribu- 
tions of these effects on the uncertainties in the IBM, we allude 
a Monte-Carlo simulation. 

The fitting algorithm in XSPEC gives, along with the best fit 
model, the covariance matrix of second derivatives calculated at 
the best-fit solution. Here, we use this covariance matrix to re- 
draw a random, new set of model BMD values on the given tem- 
perature grid (see Fig.O that is close in its statistical properties 
to the original best-fit model, in the sense that it produces spec- 
tral fits to the data that are almost as good as the best-fit model. 
In order to generate a random set of numbers corresponding to 
a correlation matrix R, one can decompose R into R - U^U. 
If X is a vector of uncorrected random values, Xc - XU is a 
set of correlated random values coiTesponding to R. The process 
of generating random BMDs is repeated numerous times, and 
for each BMD solution the coiTesponding IBMikT) is calcu- 
lated at each temperature point on the grid. Thus, we generate a 
large sample of IBM{kT) values at each temperature point, from 
which we adopt the central 68% as the equivalent of the 1-cr un- 
certainty interval that is, then, plotted in Figures [3] and I?] For 
each generated random BMD we also calculate the emission- 
measure weighted mean temperature < kT >. The central 68% 
of the < kT > distribution is adopted as the 1-cr uncertainty in- 
terval on < kT >, plotted in Figures|4]and[8] 
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